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ABSTRACT 

We study the IR emission from flared disks with and without additional optically thin halos. 
Flux calculations of a flared disk in vacuum can be considered a special case of the more 
general family of models in which the disk is imbedded in an optically thin halo. In the 
absence of such halo, flux measurements can never rule out its existence because the disk 
flaring surface defines a mathematically equivalent halo that produces the exact same flux at 
all IR wavelengths. When a flared disk with height H at its outer radius R is imbedded in a 
halo whose optical depth at visual wavelengths is Thaio, the system IR flux is dominated by 
the halo whenever rhaio > jH/R. Even when its optical depth is much smaller, the halo can 
still have a significant effect on the disk temperature profile. Imaging is the only way to rule 
out the existence of a potential halo, and we identify a decisive test that extracts a signature 
unique to flared disks from imaging observations. 
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1 INTRODUCTION 

Modeling the IR radiation of pre-main-sequence stars has tradi- 
tionally involved "classic" geometrically-thin optically-thick disks. 
This approach fails to produce many features of observed spectral 
energy distributions (SED) in both T Tauri stars (TTS) and Herbig 
Ae/Be stars (Haebes). One way out of this problem is to supplement 
the disk emission with a surrounding optically thin halo (e.g. But- 
ner, Natta & Evans 1994; Miroshnichenko et al 1999, MIVE here- 
after). A simpler alternative was proposed by Chiang & Goldreich 
(1997, CG hereafter): The surface skin of any optically thick ob- 
ject is of course optically thin. The emission from the disk surface 
layer can become significant under certain flaring conditions, and 
CG present SED fits of a number of TTS purely in terms of flared 
disks. This proposal was extended to SED modeling of Haebes by 
Chiang at al (2001) and Natta et al (2001). 

While the CG model successfully solves the SED problem 
with a simple physical explanation, recent high-resolution imaging 
observations of both TTS and Haebes reveal the presence of halos, 
ignored in the CG approach. Most striking is the case of GM Aur, 
a classical TTS whose SED modeling was presented as evidence 
of the success of the CG approach (Chiang & Goldreich 1999). 
Indeed, HST/NICMOS images by Schneider et al (2003) show a 
flared disk with radius of ~300 AU, but they also reveal the pres- 
ence of a surrounding tenuous envelope extending to a radius of 
~1000 AU. Using detailed Monte Carlo calculations, Schneider et 
al find that simplified models with flared disk without an imbed- 
ding halo fail to replicate the scattered light intensity pattern seen 
in the NICMOS images. The addition of a halo is essential for suc- 



cessful modeling of both the SED and the images. Similarly, HST 
and ground-based imaging by Stapelfeldt et al (2003) of the TTS 
HV Tau C show a nearly edge-on flared disk imbedded in a more 
extended nebulosity. From detailed model calculations Stapelfeldt 
et al find that although it is possible to fit the SED purely with a 
flared disk, the flaring is unreasonably large and such models do 
not reproduce the image adequately. They conclude that a flared 
disk alone is an inadequate model and that an additional, extended 
component to the circumstellar density distribution is needed to ex- 
plain the observations. Finally, from detailed modeling of flared 
disk SED Kikuchi, Nakamoto & Ogochi (2002) conclude that ha- 
los are necessary supplements in explaining flat-spectrum TTS. 

High resolution observations of Haebes give similar evidence 
for imbedding halos. Combining space- and ground-based observa- 
tions of HD 100546, the nearest Herbig Be star, Grady et al (2001) 
resolve the disk, which extends to 5" (~500 AU) from the star 
and displays the elongation of inclined viewing with i = 49° . But 
they also find that the disk is imbedded in a more extended halo 
(~10" radius) that is roughly circular in shape and optically thin 
(background stars are visible through it). Polomski et al (2002) per- 
formed high-resolution ground-based observations of a sample of 
Haebes whose mid-IR emission was claimed to be disk-dominated 
by Hillenbrand et al (1992). They find instead that the emission 
is not confined to an optically thick disk but originates in a more 
complex environment that includes large, extended dust envelopes. 
From analysis of ISO data, Creech-Eakman et al (2002) conclude 
that only a subset of the Haebes they observe can be described 
purely in terms of the CG model (notably, E.I. Chiang is a coau- 
thor of this study). 
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Figure 1. Model geometry and notations for a flared disk and its CG surface 
layer. Tlie unit vector 6 points at the observer's direction. 



How can the SED of the same T-Tauri star, GM Aur, be fit- 
ted successfully with a flared disk both with a halo (Schneider et al 
2003) and without one (Chiang & Goldreich 1999)? Is it at all pos- 
sible to distinguish between these two cases? What, if any, is the 
unique radiative signature of each component? Obviously, imbed- 
ding a disk in a tenuous halo with a very small optical depth is not 
going to affect its emission appreciably. At which stage then does 
the halo assume a significant role? Here we address these questions. 
In appendix A we derive general results for the emission from op- 
tically thin dust in an arbitrary geometry. We apply these results to 
the surface layers of flared disks in §2 and to halos in which the 
disks could be imbedded in 53. 



2 FLARED DISKS 

2.1 The CG layer and its equivalence to spherical halo 

Chiang and Goldreich (1997) noted that the emission from the op- 
tically thin surface layer of an optically thick disk, which has been 
neglected traditionally, can become significant under certain flaring 
conditions. The stellar radiation penetrates to an optical distance rv 
= 1 along a direction slanted to the surface by angle a (figure 0. 
The optical depth of the corresponding skin layer along the nor- 
mal to the surface n is a at visual and aqv at wavelength v, where 
Qv = cTy/crv and ctv is the dust cross-section at visual. In the case 
of a flat thin disk whose inner radius is determined by dust subli- 
mation, the grazing angle is 



a* , . 4 i?,^ 

aflat = — , where a = — -?r> 
a i-K Kb 



(1) 



here i?^, is the radius of the star, is the dust sublimation radius 
(eg. lAlt and ra is distance from the axis. Flaring is defined by the 
radial profile of the disk height H ra) or, equivalently. 



H H 
p = arctan — ~ — . 



(2) 



As is evident from figures a = 7 — /3 where tan 7 = dH/dra_, 
therefore the grazing angle of a flared disk is 

dl3 



da 



(3) 



The CG surface layer serves as an effective optically thin disk 
atop the underlying optically thick disk core, and its flux is ob- 
tained from the volume integration listed in equation IA7I In the 



case of face-on orientation, the optical depth of the CG layer obeys 
o'vj n^dz — q,^a, therefore its flux at distance D is 



2-KRt 



/ Bi,{T) ctada. 



(4) 



Since the temperature profile of optically thin dust depends only 
on distance from the radiation source, the geometry dependence of 
this expression enters only from the radial variation of the graz- 
ing angle a (reflecting the dust column variation). But other ge- 
ometries can produce an identical expression. For example, in the 
case of spherical geometry the flux is controlled by the dimension- 
less density profile r) (equation lASt and the radial optical depth is 
Tu = <Jij^ Hddr. Denoting the optical depth at visual wavelengths 
by TV, the flux is 



sph, i/ 



47tRI 

D2 



g^rv / B^(T)Tiy^dy, 



(5) 



where y = r/Rs. Since y and a enter only as integration variables 
in the last two integrals, the two expressions are mathematically 
identical if 



a(y) 

■q oc — and rv 



da 



(6) 



A minor point is that T in eq. |4] is strictly a function of y = 
ayT+"^ rather than a (since temperature is controlled by dis- 
tance from the star); this slight difference can be ignored because 
/3 ^ 1 everywhere in the disk. 

Scattering can be treated similarly and produces the same re- 
sult since the only difference is that Bv (T) is replaced by = 
L^/A-Kr"^ , another function of y (see equation lASt . Therefore, there 
is a complete equivalence between spherical halos and the CG sur- 
face layers of flared disks; either case defines a model with the other 
geometry and the exact same flux. Specifically, equations Q and 
|6| show that a thin flat disk is equivalent to a spherical halo with 
r; oc l/y^ and rv = ^Q*- In the case of a flared disk, equations |3| 
and|6|show that the CG layer will produce precisely the same flux 
as a spherical halo with 



77 cx 



d/3 
dy' 



rv = i[/3(i?d)-/3(i?s)]; 



(V) 



in particular, the equivalent halo of a disk with flaring angle /3 cx 
1/a^ has density distribution 77 oc . And conversely — given 

a spherical halo, the flared disk with the same outer radius and 



/?(a) =/3(l) + 2rv / n{y)dy 



(8) 



will produce the exact same flux from its CG-surface layer. 

Each disk defines a mathematically equivalent halo. Although 
this equivalence was derived only for face-on viewing of the disk, 
it carries to most inclination angles since the flux from optically 
thin dust involves a volume integration (equation IA7> and the ob- 
served fraction of the disk surface layer remains largely intact as 
long as internal occultation is not too significant. Similarly, the gen- 
eral analysis presented in appendix A shows that the spherical ide- 
alization is not essential for the halo geometry. The dust distribu- 
tion can be flattened and even distorted into irregular shape before 
severely affecting the emerging flux. 

These results resolve the paradox of successful SED fits for the 
same star, GM Aur, with a flared disk both with a halo (Schneider 
et al 2003) and without one (Chiang & Goldreich 1999). The halo 
contribution to the flux can be absorbed into the disk component by 
readjusting the flaring law, enabling a successful SED fit without a 
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Figure 2. Points on the surface of a flared disk at equal distance from tlie 
star lie on a circle centered on the disk axis. The circle retains its shape in 
pole-on viewing but is deformed into an off-center ellipse (eguation ll II in 
viewing from inclination angle i. 



halo even though it is directly visible in imaging observations. The 
problems Stapelfeldt et al (2003) report with disk-only models of 
HV Tau C find a similar explanation. The halo contribution to the 
SED can be shifted to the disk, with the halo optical depth added 
to the flaring angle (equation 0. This leads to unreasonably large 
flaring, as Stapelfeldt et al find. The same problem was encountered 
by Kikuchi et al (2002) in modeling flat-spectrum TTS. The halo 
optical depth in these sources are sufficiently high that they cannot 
be substituted by disks with realistic flaring. 

In addition to explaining the shortcomings of disk-only mod- 
els in these specific cases, the equivalence between halos and CG 
layers has two important broad consequences: 

(i) When the disk is imbedded in a halo that radiates more than 
its CG layer, the halo becomes the dominant component of the IR 
flux. This happens when the halo contains more dust than the halo- 
equivalent of the disk. From equations |2| and |7| the IR radiation 
from the system is dominated by the halo contribution whenever 

nialo > 7^ (y) 

4 Hd 

where Thaio is the optical depth across the halo at visual wave- 
lengths (2tv for spherical halos). 

(ii) It is impossible to distinguish the CG layer of a flared disk 
from a halo with flux measurements. Only imaging can produce an 
unambiguous signature of the CG layer. 



2.2 Images 

The brightness contours of face-on flared disks are concentric cir- 
cles centered on the star. Inclined viewing changes the contours 
substantially. Consider the intensity of radiation scattered from the 
CG surface layer. It obeys 



Joe 



r(6) 



(10) 



where r is the distance to the star and r(o) the optical depth toward 
the observer at the scattering point (equation |A2j. Both factors in- 
troduce distinct image asymmetries. 



The fundamental reason for image distortion by inclination is 
that the same projected distance from the star corresponds to widely 
different locations on the surface of the disk. On that surface, con- 
tours of equal distance from the star are circles of radius r^. When 
viewed face-on from distance D, each contour appears as a con- 
centric circle of radius = r^,/ D, as seen in the top image in 
figure|2| At inclination viewing angle i to the disk axis, the contour 
is no longer circular. Absent flaring, the contour becomes an ellipse 
centered on the star with major axis 2^a and minor axis 2Sa cos i, 
aligned with the projection of the disk axis on the plane of the sky. 
Flaring raises the contour to height H = tan /3 above the equa- 
torial plane (figure 0, and the star is shifted toward the observer 
along the minor axis by 6^ tan (3 sin i. A point on the contour at 
position angle (j) from the near side of the minor axis is observed at 
displacement 6 — 9ag{<j)) from the star, where 

g{<f)) — [(tan /3 sin i — cos i cos (f>)'^ + sin^ 0] ^''^ . (11) 

These contours are shown in the bottom image of figure|2| The off- 
center position of the star on the minor axis creates an asymmetry 
such that the far and near portions of this axis obey 6'far/Sncar ~ 
cos(i — 13)/ cos(i + 13). Because (3 increases with ^a, this asymme- 
try increases with distance from the star. 

At observed displacement [6, (f>) from the star, a point on the 
surface of the disk is located at r ~ r-a = D9/g{(j)). At that point 
the optical depth of the CG layer toward the observer is r(6) = 
q,^a/o{<f)), where 

o{<f)) — n ■ 6 — cos i cos 7 — sin i sin 7 cos 0. (12) 



Therefore the scattering image obeys 



5(0) 



0(0)- 



(13) 



In this expression, both a and 7 are determined at the location 
Ta = D6/g(<j)) on the disk surface. A power-law grazing angle 
Q oc produces the image I{6,(f)) oc [g{<j>)/6]'^'^'' /o(4>). This 
expression and the resulting brightness contours explain the results 
of Monte Carlo scattering calculations for flared disks (Whitney & 
Hartmann 1992, Wood et al 1998). Figure |3| shows the scattering 
images at three viewing angles of a flared disk with the parameters 
suggested for AB Aur by Dullemond et al (2001). 

The images produced at emission wavelengths are handled in 
complete analogy. The only change is the replacement of r^^ by 
the temperature T, i.e., another function of r, modifying the de- 
pendence of brightness on g[(j>)/6. 



2.2.1 Image asymmetry 

Brightness contours not subject to rim occultation are ellipses with 
eccentricity e = cos i that directly determines the inclination angle 
irrespective of the flaring profile. The images shown in figure |3| 
possess an additional deviation from circular symmetry, unique to 
flaring and conveniently measured by the brightness at diametric 
locations across an axis through the star 

i{e,4> + ^)^i{e,^) 



A{e,d)) 



(14) 



i{e,4> + i,) + i{e,^y 

This asymmetry parameter vanishes for flat disks at all inclination 
angles, and for pole-on and edge-on viewing irrespective of the flar- 
ing. However, at intermediate inclination angles, flaring introduces 
substantial asymmetry, as is evident from figure|3| 

Non-vanishing A is the hallmark of inclined flared disks be- 
cause it measures the displacement of the isophote centers from the 
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Figure 3. Viewing at inclination angles i = 25°, 45° and 65° of radiation 
scattered off the surface of a flared disk with 3" radius and constant grazing 
angle a = 4.9°. For each i, the left panel shows the image with brightness 
contours, the right panel the asymmetry factor A (eguation lMI along three 
azimuthal directions, shown in the top left panel, whose angles are des- 
ignated from the observer's direction. The asymmetry parameter vanishes 
along the major axis (</> = 90°) at all inclination angles and along any axis at 
i = 0° and 90°. For flat disks, A is identically zero at all inclination angles. 



peak brightness position. Its systematic variation with azimuthal 
angle easily distinguishes it from deviations from the perfect ge- 
ometry of idealized models or noise in the data. Each flaring profile 
produces its own characteristic signature A. For example, it is easy 
to show that the constant grazing angle a used in figure |3| gives 
A ~ tan (3 tan i along the minor axis. Therefore, measuring A de- 
termines the flaring profile once the inclination is determined from 
the eccentricity of the brightness contours. 



3 HALO-IMBEDDED-DISKS 

The results of the previous section show that flux calculations of 
a flared disk in vacuum can be considered a special case of the 
more general family of models in which the disk is imbedded in 
an optically thin halo. In the CG case the "halo" is the disk surface 
layer, fully determined from the flaring geometry. This layer pro- 
vides the same IR emission and heating of the underlying optically 
thick core as its equivalent halo. Therefore, a study of the general 
halo-imbedded-disk problem contains every possible CG model of 
flared disks while also covering all cases in which the disk is indeed 
imbedded in a dusty halo whose optical depth exceeds the bound in 
equation|9| 

Consider a star surrounded by a geometrically thin passive 
disk and a dusty halo (figure|4j. We study the case of a flat disk and 
a spherical halo, simplifications that enable us to derive analytical 




Figure 4. Geometry of the halo-imbedded-disk model: a flat geometrically- 
thin optically-thick disk extends from the stellar surface to radius . An 
optically thin spherical halo extends from the dust sublimation radius Rs to 
-Rh . The small pillbox at the disk surface serves as a Gaussian surface for 
flux conservation. 



results and broad conclusions that offer important insight. These 
simplifications do not cause any serious limitations. As noted pre- 
viously, the halo geometry can be distorted considerably without 
much impact on the outcome. And because only the optically thick 
core of the disk enters into considerations, its only relevant prop- 
erty is its temperature distribution, the surface shape is immaterial. 
Indeed, Wolf et al (2003) present detailed model calculation of a 
flared disk imbedded in an envelope pinched around the equatorial 
plane, and their results fully conform to our conclusions. 

The halo extends from the inner radius Rs to some outer radius 
-Rh = Y Rs - Thanks to scaling, instead of these radii we can specify 
the dust temperature on each boundary (see appendix A). The halo 
is fully characterized by its density profile r;(j/) (eg. lASt and optical 
depth rv; only rv < 1 is relevant in TTS and Haebes since the star 
is always visible. Because of its potentially large optical depth, the 
disk can extend inside the dust-free cavity where its optical depth 
comes from the gaseous component. The geometrically-thin disk 
assumption implies that the disk temperature varies only with ra- 
dius, vertical temperature structure is ignored. This temperature is 
calculated from radiative flux conservation through the Gaussian 
surface in the shape of a small pillbox shown in figure |4| Denote 
by Ti the radiative flux entering the pillbox from above, including 
both the stellar and diffuse components. Then 

jB,,{T)E3{T°)di' + Jfihe'^^^^dfldu = H, (15) 

27r 

where E3 is the 3rd exponential integral, is the disk vertical 
optical depth and fj, is the cosine of the angle from the disk nor- 
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mal. The first two term on the left are the disk contribution to the 
upward flux, the third is the contribution of local intensity trans- 
mitted upward through the disk. When the disk is optically thick at 
frequencies around the peak of B„{T) (r° ^ 1 for ~ kT/h), 
the second and third terms can be neglected, leading to the stan- 
dard expression for disk temperature (e.g. Friedjung 1985, Kenyon 
& Hartmann 1987). We assume this to be the case everywhere in 
the disk, an assumption that we check for self-consistency in all our 
model calculations. 

With its temperature derived, the disk emission is calculated 
from By (r)[l — exp(— r,^ /m)]- Because the disk is optically thick 
around the Planckian peak at all radii, the self-absorption factor 
can be neglected in the calculation of the disk overall bolometric 
flux. Then Fdisk{D, i), the disk flux observed at distance D and 
inclination i, is proportional to cos i, reflecting the variation of pro- 
jected area. This proportionality remains largely unaffected by the 
envelope attenuation because the short wavelengths, the main con- 
tributors to the bolometric flux, emanate from close to the star so 
that their pathlength is approximately isotropic. Denote by I/disk 
the disk contribution to the overall luminosity L and by Lsph the 
contribution of the (halo -I- attenuated stellar) spherical component. 
The corresponding flux components are then 

-^sph 



1000 



r-1 / i~\ ■\ -^disk 

-Fdisk(-D,«) = -^i^ cos I, 
and the overall flux is 



sph 



47rD2 



(16) 



F{D, i) = Fdisk(-D, 0) cos I + F,pi,{D) 



1 + 2x cos i 
1 + 2: 



where x — -Ldisk/isph. The standard "bare" disk has Ldisk = \L 
(Kenyon & Hartmann 1987), therefore in this case x — ^. Larger 
fractions can occur when the disk is imbedded in a halo because of 
the heating effect of the diffuse radiation, discussed below. 

We performed detailed model calculations with the code 
DUSTY (Ivezic, Nenkova & Elitzur 1999) which takes into account 
the energy exchange between the star, halo and disk, including dust 
scattering, absorption and emission. Because its optical depth is 
typically tv < 1, the halo is transparent to the disk emission in 
all the models we consider and we neglect the disk effect on the 
halo. In all the calculations, the spectral shapes of the grain ab- 
sorption and scattering coefficients are those of standard interstellar 
mix, the sublimation temperature Ts = 1500 K. The spectral shape 
of the stellar radiation is taken from the appropriate Kurucz (1994) 
model atmosphere. 



3.1 Temperature Profiles 

The heating rate of a thin flat disk by the stellar radiation at r 2> i?* 
is 



n. 



2fs 1 
Stt Rs ' 



(17) 



where Fb is the stellar flux at Rs and a = r^/Ra, with ra distance 
from the axis (Friedjung 1985). This result reflects the decline 
of the stellar solid angle and the 1/a variation of the grazing angle, 
yielding disk temperature variation T oc a"'^''*. Natta (1993) noted 
that imbedding the disk in a dusty halo can significantly affect its 
temperature even at small halo optical depths (see also D'Alessio, 
Calvet & Hartmann 1997 for the effect of optically thick halos). 
With a simple model for scattering at a single wavelength she found 
that the disk temperature law becomes T oc q^^^+p)/* if the halo 
density profile is oc y^^. 
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Figure 5. Top: Temperature profiles of a disk when heated only by a central 
star with = 10,000 K (full line), and when imbedded in a spherical dusty 
halo with tv = 0. 1 or 1 , as marked. The halo starts at dust sublimation Ts 
= 1,500 K and its density profile is oc r"^ The temperature profile of the 
halo is also shown in each case. Bottom: The fractional contributions of the 
halo and (attenuated) stellar components to heating of the disk 



Our calculations confirm this important point. Figure|5|shows 
the temperature profile for a disk around a T* = 10,000 K star when 
"bare" and when imbedded in a spherical halo with i] oc and 
TV = 0.1 and 1^. Even though a halo with rv as small as 0.1 is 
almost transparent to the stellar radiation, it still causes a large rise 
in disk temperature. As is evident from the bottom panel, the halo 
contribution to the disk heating overtakes the stellar contribution 
inside the dust-free cavity and dominates completely once the dust 
is entered. Wolf et al (2003) present a similar figure for their model. 

A dusty envelope with ry = 0. 1 intercepts only about 10% of 
the stellar luminosity while the disk intercepts 25% of that lumi- 
nosity. So how can the halo dominate the disk heating? The reason 
is that direct heating of the disk by the star occurs predominantly 
at small radii. The disk absorbs more than 90% of its full stellar 
allotment within lOi?* while its entire remaining area, even though 
so much larger, absorbs only 0.0251/. From equation IA II the halo 
starts at Rs ~ 1007?*, where Tit has already declined to ~ 10~® 
of its value at the stellar surface. In contrast, the halo emission is 
isotropic, therefore half of the radiation it intercepts is re-radiated 



^ The addition of a halo can only raise the disk temperature, yet figure 151 
shows that our calculated profile for Ty= 1 dips slightly below that of the 
bare disk at a < 0.2. This happens because we neglect the disk emis- 
sion in the calculation of the halo temperature. The eiTor introduced by this 
approximation is of order 2% when rv= 1, and even less at smaller ry. 
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Figure 6. Sample SEDs for halo-imbedded-disks around stars with T* = 10,000 K. The halo starts at Rs where Ts = 1,500 K and extends to l,OOOHs, with 
density profile and optical depth rv as indicated. The disk starts at the stellar surface and extends to the radius set by the edge temperature T^j = 25 K 
(see Table 1). Each SED is normalized with the bolometric flux Pboi = i) of the appropriate viewing angle (see eQuation ll6l . Top row: The SED for 

pole-on viewing. Full hues denote the overall flux, dotted Unes the spherical (halo + attenuated stellar) component and dashed lines the disk component. The 
thick dot-dashed line is the flux from a face-on "bare" disk in the absence of an imbedding halo. It is omitted in the other rows. Mid row: Same models viewed 
at inclination angle i = 85°. Bottom row: Variation of the overall SED with viewing angle i. Results for i < 70° are barely distinguishable from i = 0°. 



toward the disk, greatly exceeding the direct stellar contribution. 
Further insight can be gained from the approximate solution pre- 
sented in Ivezic & Elitzur (1997; IE hereafter) for radiative transfer 
in spherical symmetry. From equations 20 and B4 of IE it follows 
that disk heating by a halo with rv < 1 and rj oc is roughly 



rth = — — rv X 

8 p + 1 



1 



1/a 



o < 1 
a > 1 



(18) 



And since Rs ~ lOOii*, the halo dominates the heating at a = 1 for 
rv as small as 0.02 when p = 2. As rv increases, the halo domi- 
nance of the heating moves inside the cavity. There the halo heating 
remains approximately constant while the stellar heating varies as 



therefore stellar heating dominates only at a < (60rv) 



-1/3 



at larger distances the halo takes over. This explains the results pre- 
sented in the lower panel of figure |5| as well as figure 9 in Wolf et 
al (2003). 



when p > 1; when p = 1 the factor (p — l)/(p + 1) is replaced 
by 1/(2 In y). This yields T oc a~^^+^^'^, corroborating Natta's 
result and extending its validity beyond the single-wavelength scat- 
tering approximation she employed. The temperature profile is sim- 
ilar to that of a bare disk when p — 2 but more moderate when 
p < 2. More importantly, 



97r p - 1 i?s 
T6p+ 1 % 



■ TV- 



(19) 



The figure also shows the temperature profile of the halo. This 
profile is largely independent of rv, varying roughly as 
when the long-wavelength spectral shape of the dust absorption co- 
efficient is qi, cc v"" . The important property evident in the figure is 
that the disk is much cooler than the envelope at all radii at which 
both exist and can also contain cooler material in spite of being 
much smaller, with far reaching consequences for the system radi- 
ation. 
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Tv= 0.1 Ty= 1 



X 


Yd 


X 




0.43 


190 


1.33 


400 


0.47 


85 


1.56 


135 



Table 1. Derived parameters for the models whose SEDs are presented in 
figure|S| The luminosity ratio of the components is x = idisk/^sph ^nd 
the disk radius is = y^/Jg, set from the requirement = 25 K. A "bare" 
disk (ry = 0) has a; = ^ and = 18. 

3.2 SED 

From eguation fT^l the fractional contribution of the disk to the over- 
all bolometric flux is 



-Fdisk + Fsph 1 + 2x cos i ■ 

Face-on orientation gives the maximal p = 2x/{l + 2x) and the 
standard "bare" disk, with x = \, has n C t- Introduce the nor- 
malized SED = F^/ \ F^dv, with similar, separate definitions 
for the disk and spherical components of the flux. Then 

fv = p/i/,disk + (1 — P).fv.sph (21) 

Since the disk flux obeys Fi/,disk(i) ~ -Fj^,disk(0) x cosz for the 
range of parameters considered here, /^,disk is independent of the 
viewing angle i, and the entire i-dependence of the SED comes 
from the mixing factor p. 

Figure |6| shows sample SEDs for some representative mod- 
els. The halos extend from Rs to l,000i?s, with density profiles 
and overall optical depths as indicated. The behavior of SEDs for 
spherical shells was discussed in IE and since the halo emission is 
unaffected by the imbedded disk, the SEDs plotted in dashed lines 
need no further discussion. The disk, on the other hand, is strongly 
affected by the halo as is evident from contrasting each disk SED, 
plotted in long-dashed line, with what it would have been in the 
absence of a halo (dot-dashed line). The two curves are identical 
within the first bump around I ^m, caused by the stellar heating. 
In the absence of a halo, the disk SED drops from that peak as 
XFx oc A"*''"^. However, halo heating of the outer regions of the 
disk generates the second broad bump of disk radiation, which is 
almost two orders of magnitude higher than the "bare" disk emis- 
sion. 

The halo heating effect is also evident from other disk proper- 
ties. The disks in these models start at the stellar surface and extend 
to a radius where the temperature is Td = 25 K. In the absence 
of a halo, this temperature would be reached at Vd = Rd/Rs = 18. 
As Table 1 shows, heating by even a tenuous halo with rv = 0. 1 
pushes this radius out by almost factor 5 for the steep density pro- 
file and another factor of 2 for the flatter r^^ profile which 
spreads the heating further away from the star. The impact of halo 
heating increases with rv, pushing Rd further out still. Similarly, 
the disk luminosity, xL/(l + x), is only 0.251/ in the absence of 
a halo but increases as the halo directs more radiation toward it to 
the point that it becomes 0.61/ when rv = I. 

Although it is more compact, the disk can become the stronger 
emitter at long wavelengths so that the overall SED is dominated 
by the halo at IR wavelengths and by the disk at sub-mm and mm 
wavelengths. This role reversal affects also the wavelength behav- 
ior of images. Figure shows a series of images at various wave- 
lengths for a sample model. At IR wavelengths the image is dom- 
inated by the halo, displaying the size variation discussed in ap- 
pendix A. The 0.6pm image is dominated by scattering, the A ^ 




Figure 7. Images at various wavelengths of a halo-imbedded-disk around 
a star with = 10,000 K. The halo starts at Rs where % = 1,500 K and 
extends to l,000_Rs, with density profile r] oc y^^ + O.Oby^^ and rv = 
0.6. The disk starts at the stellar surface and extends to radius Rd set by 
= 25 K. The angular scale corresponds to bolometric flux -Fbol = 10"^" W 
m~^. The viewing inclination angle is j = 76°. 

5pm images reflect dust emission, leading to size increase with 
wavelength. The disk emerges at llpm and dominates the A ^ 
100/im images. The finite beam size and dynamic range of any 
given telescope could result in an apparent size decrease between 
10 pm and 100 ^m in this case. Such an effect has indeed been dis- 
covered in the Haebes MWC 137, whose observed size decreases 
between 50 pm and 100 pm (DiFrancesco et al 1994, 1998). A 
switch from envelope to disk domination provides a simple ex- 
planation for this puzzling behavior. No single dust configuration 
can explain such a decrease, a conclusion reached already in MIVE 
and further affirmed by the results of appendix A. DiFrancesco et 
al (1998) suggest that this behavior might reflect multiple, rather 
than singular, sources of heating but the results of appendix A show 
the inadequacy of this conjecture. The region heated by any single 
source displays an increase of observed size with wavelength, and 
the superposition of multiple heating sources preserves this behav- 
ior. The opposite trend is possible only when the density distribu- 
tion contains two distinct components, one optically thick, cool and 
compact, the other optically thin, warmer and more extended. 

A similar effect was detected also in the dust-shrouded main- 
sequence star Vega. Van der Bliek, Prusti & Waters (1994) find that 
its 60 pm size is 35" ± 5", yet 850 pm imaging by Holland et al. 
(1998) produced a size of only 24x21" ± 3". So the dust distribu- 
tion around Vega, too, could contain both spherical and disk com- 
ponents with the switch from halo- to disk-dominance occurring 
somewhere between 60 and 850 pm. Indeed, imaging at 1.3mm by 
Willner et al (2002) reveal the presence of the disk. 

The calculation of the disk emission contains two free param- 
eters, involving the temperature and optical depth at the disk edge. 
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Figure 8. The effect on the SED of varying the temperature (T^j) and 
350/xm optical depth (T350) of the disk outer edge. The displayed model 
has rv = 0.5 and rj oc r~^. The viewing angle is 85°. Dashed lines cor- 
respond to r35o = 1, full lines to T350 sufficiently large that the disk edge 
remains optically thick at all displayed wavelengths. 



The disk outer radius i?d determines its lowest temperature Td 
and a corresponding Planck-peak wavelength. While shorter wave- 
lengths are emitted from a range of disk radii, all longer wave- 
lengths originate from the edge of the disk. The resulting effect 
can be seen in figure|8| which shows two representative values of 
Td. By the model assumptions the disk must be optically thick at 
the local Planck-peak wavelength everywhere, ~ 350 at Td 
= 25 K. As long as the disk edge remains optically thick also at 
longer wavelengths, the emission follows the Rayleigh-Jeans pro- 
file /^.disk oc . Once the disk edge becomes optically thin, the 
SED switches to the steeper decline /^,disk oc v^Ov at longer wave- 
lengths. The break in the disk SED is controlled by the optical depth 
of the disk edge, which we specify at 350 /xm. Figure|8|shows also 
the effects of varying raso from its smallest value rgso ~ 1 to a 
value sufficiently large that the edge is optically thick at all the dis- 
played wavelengths. 



3.2. 1 The disk inner radius 

The bottom panel of each model in figure |6| shows the variation of 
the overall SED with viewing angle i. The entire variation comes 
from the mixing factor p (see equations 1201 and 1211 . Since the pa- 
rameters X and i enter only in the product x cos i but not separately, 
the SED is subject to a degeneracy: systems viewed at different 
inclination angles will have the same SED if they have the same 
X cos i in addition to all other properties. Because of the rapid de- 
cline with distance of the radiation absorbed from the star (cf eq. 
I17> . the disk luminosity, i.e., x, has a steep dependence on its in- 




R,/R. 

Figure 9. Contours of fixed mixing coefficient p (eguation llOL as marked: 
The SEDs are the same as those presented in figure|^for i = 70° and 85° 
when the disk inclination i and its inner radius /?in vary together along each 
of the plotted curves. The contours are virtually the same for the and 
halo density profiles. 



ner radius i?in; moving the disk inner edge from to only IRi, 
removes 56% of the stellar luminosity intercepted by the disk, 3-R* 
results in a 72% removal. Such central holes reduce x but do not 
impact any other relevant property because they remove only the 
hottest disk material whose contribution to the overall flux is negli- 
gible in comparison with the stellar component. 

Figure|9|plots contours in the i-i?in plane of constant mixing 
factor p. It shows, for example, that the SEDs presented in figure 
|6|for i = 70° would be the same for systems viewed at z = 35° 
if the disk inner radius is increased from li?* to 2.2_R* in the rv 
= 0. 1 case and 6-R* for rv = 1 . Although the sizes of these holes 
cannot be determined from SED modeling of single stars, from sta- 
tistical arguments MIVE conclude that their existence is essential 
to produce a plausible distribution of inclination angles. 



3.2.2 Flared disks without halos 

Thanks to the CG layer-halo equivalence, each model presented 
here describes also a disk with no halo and with the flaring angle 
defined by eq. |8|from 77 and rv. In particular, the SED presented 
for rj oc halos describe also flared disks without halo and 

with a oc 1/a, those with rj (x 1/y cover flared disks with constant 
grazing angle a (/? ~ In a). The equivalence strictly holds only for 
disks and halos of the same size. However, the models presented 
would be little affected if each halo was truncated at i?d because 
that would only remove halo emission at long wavelengths where 
the SED is dominated by the disk anyhow. Since the flaring angle of 
each equivalent disk reaches 2rv at its outer edge, halos whose op- 
tical depths require excessive flaring cannot be realized with disks 
only. This was the problem recognized by Stapelfeldt et al (2003) 
and Kikuchi et al (2002). 

The equivalence holds only for the SED. High-resolution ob- 
servations would produce widely different images for each halo and 
its equivalent flared disk, except when the latter is viewed face-on. 
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4 DISCUSSION 

The results of il2. II show that every disk imbedded in a halo with 
TV < I can be replaced by a flared disk without a halo and with an 
identical flux. This mathematical degeneracy explains the success 
of SED modeling with disk alone a system like GM Aur (Chiang 
& Goldreich 1999) even though the halo was subsequently discov- 
ered in imaging observations (Schneider et al 2003). It also explains 
why such modeling runs into difficulties and requires excessive flar- 
ing when the halos have larger rv, as is the case in HV Tau C 
(Stapelfeldt et al 2003) and flat-spectrum TTS in general (Kikuchi 
et al 2002). 

Disks generally do not exist in pure vacuum. Equation |9| de- 
fines the circumstances under which the surrounding dust becomes 
the dominant component of the IR flux. Even at smaller rv, when 
not dominating the overall flux, the halo can still dominate the disk 
heating and make a strong impact on its temperature profile ( i)3.1> . 
Ignoring the surrounding material can produce misleading results 
regarding the disk properties, such as its flaring profile. 

In spite of the attractiveness of the flared disk as a simple, 
physical model without any additional components, imaging ob- 
servations give irrefutable evidence for the existence of extended 
halos in pre main-sequence stars. The origin of these halos has not 
been studied yet. Stapelfeldt et al (2003) suggest that a replenish- 
ment process, either continued infall from the surrounding ISM or a 
dusty outflow from the source itself, is operating. It is noteworthy in 
this regard that accretion with the small rates of ~ 10~* Mq yr^^ 
has been deduced from UV spectra of both Haebes (Grady et al. 
1996) and TTS (Valenti, Basri & Johns 1993; Gullbring et al 1998; 
Lamzin, Stempels & Piskunov 2001), and is consistent with halos 
that have tv ~ 0. 1 (MIVE). These low rates cannot correspond 
to the main accretion buildup of the star but rather a much later 
phase, involving small, residual accretion from the environment. 
The corresponding accretion luminosities are only ~ 0.1 Lq, justi- 
fying their neglect in our calculations. 

In addition to the CG layer-halo equivalence, our results re- 
veal numerous degeneracies that underscore the severe limitations 
of attempts to determine the dust morphology from SED analysis 
without imaging observations. The SED of a halo-imbedded-disk 
remains the same when the viewing angle and the size of the disk 
central hole vary together as shown in figure|9| From the results of 
appendix A, the SED of a spherical shell with power-law density 
profile displays a dependence on p only when p < 2.6, and 
then only in the spectral region A ^ Aout . All other regions of p and 
A produce the same universal behavior Fy oc u^Ou ( i|A4> . The dust 
optical properties introduce additional degeneracies. The results of 
i|A4l show that the frequency dependence of Gi, and the radial de- 
pendence of the density profile can be interchanged on occasion 
without affecting the SED. The fundamental reason for all these 
degeneracies is that the flux from an optically thin source involves 
a volume integration (eq. IA7t that tends to remove much of the 
dependence on the underlying morphology. Although the specific 
degeneracies we uncovered involve spherical geometry, the general 
analysis in appendix A shows that the spherical idealization is not 
essential. The dust distribution can be flattened and even distorted 
into irregular shape before severely affecting the results. 

These degeneracies make it impossible to determine the ge- 
ometry from a fit to the SED alone without additional input. Only 
imaging can trace the actual density distribution, and scattering pro- 
vides a more faithful presentation because, unlike emission which 
involves also the dust temperature, it involves only the density dis- 
tribution. Reliance on SED modeling alone can produce misleading 
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results, as was the case for the parameters of the flared disk in GM 
Aur. 
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APPENDIX A: OPTICALLY THIN DUST 

Consider a cloud heated from inside by a star of radius i?* and 
effective temperature T^, (figure lAll . The star clears out a dust-free 
cavity of radius Rs, determined by dust sublimation T{Rs) = Ta. 
When the dust is optically thin, the cavity radius can be found from 
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out 



Figure Al. A star imbedded in a cloud clears out a dust-free cavity of radius 
-Rs, corresponding to dust sublimation temperature Ts (eguation lAU . The 
dust temperature declines with radial distance toward its surface value Tout ■ 
The intensity at impact parameter b is obtained from integration along the 
indicated path toward the observer. 



where ^{T) is the Planck average at temperature T of the absorp- 
tion cross section (IE). With standard interstellar dust and Ts = 
1500 K, the cavity radius obeys Rs/R* — 100 at a typical Haebes 
temperature T* = 10,000 K. In ITS, on the other hand, the dust is 
much closer to the star: 7?s/i?* is only 15 at = 5,000 K and as 
small as 3 at = 3,000 K. The intensity at frequency u and impact 
parameter b (fig. lAH is 



Lib) 



(1 — ■(JJv)Bv + ■(^UvJv] rid dz, 



(A2) 



assuming isotropic scattering. Here na is the dust density, z is 
distance along the path to the observer, vj^ is the albedo at fre- 
quency u and Jv = J Iudfl/4:-K. This expression neglects self- 
absorption by the dust; the error in this approximation is of order 
1 — exp (— J (j^Uddz^. 



Al Scattering Wavelengths 

Since the dust temperature cannot exceed the sublimation temper- 
ature Ts, there is no dust emission at A < 4/imx(1000K)/rs, only 
scattering. Diffuse radiation and attenuation between the star and 
the scattering point can be neglected since our discussion is cen- 
tered on optically thin dust. Then the only source of scattering is 
the stellar radiation with energy density = Lv/A-kt'^ , where 
Ly is the stellar luminosity at frequency u. From eQuation lA2l the 
scattered brightness is 



, > f Uddz 

47r / 



(A3) 



where Q = hjD and D is the distance to the observer. Since 
the frequency- and geometry-dependence separate out, all scatter- 
ing wavelengths share a common image. Only the brightness level 
varies with v, and because of the wavelength decline of vj^Ov the 
observed size generally decreases with wavelength when traced to 
the same brightness level. In any geometry the scattering image al- 
ways traces directly the variation of column density along the line 
of sight, the dust temperature profile is irrelevant. 



A2 Emission Wavelengths 

At wavelengths longer than '~3/im, zui, < 10^^ and scattering can 
be neglected. The Planckian enters in equation I A2I as a function of 
T at fixed v, which can be well approximated by its Rayleigh- Jeans 
limit atT ^ = Q.hQhv/k and a sharp cutoff at T^. With this 
approximation the integration is limited to locations along the path 
where T ^ T,y, regions with T < T^, are too cold to emit apprecia- 
bly at frequency v. Excluding highly patchy geometries, the highest 
temperature on the path occurs at r = & (i.e., z = 0), the closet dis- 
tance to the star, and only paths with T{b) > contribute to the 
brightness. As z increases in either direction, T decreases. The in- 
tegration is truncated either because the temperature becomes too 
low, in which case the emission is temperature bounded, or be- 
cause the edge of the source is reached and the emission is matter 
bounded. Denote the resulting integration limits Zi (i = I, 2), then 



Z2 



kT rid dz. 



(A4) 



In the matter bounded case Zi is the edge of the source, the integral 
is independent of f and the frequency dependence of the inten- 
sity follows u^a^. In the case of temperature bounded emission the 
integration limits introduce additional i^-dependencies that modify 
this behavior However, the integration can be extended to oo when- 
ever (I) Zi ^ b and (2) the product UdT of dust density and tem- 
perature declines along the path faster than 1/z. Therefore, when 
these two conditions are met, the frequency variation of optically 
thin emission is Iv oc v^av even when it is temperature bounded. 
Independent of geometry, all frequencies that obey these conditions 
produce a common image, similar to the scattering case; only the 
scale of brightness varies with v. This result makes it possible to 
determine the wavelength dependence of the dust cross section di- 
rectly from imaging observations. 

Similar to the variation along the line of sight, when b in- 
creases (moving away from the star) the emission again is truncated 
by either the matter or temperature distribution. Denote by Tout the 
temperature at the source outer edge. The corresponding emission 
cutoff wavelength is 



Ao 



lOO/imx 



40K 

Tout 



(A5) 



When A > Aout , the dust is sufficiently warm everywhere that the 
emission is truncated only by the matter distribution. The observed 
size is then Q, the angular displacement of the source edge from 
the star, the same for all wavelengths. However, when A < Aout 
the brightness is truncated when T{b) ^ before the edge of 
the source is reached, resulting in a wavelength-dependent angular 
size Ox < O. The observed size of optically thin emission increases 
with wavelength so long as A < Aout, the opposite of the trend at 
scattering wavelengths. 

The frequency variation of the dust cross section is well de- 
scribed by Ov oc with n = 1-2. Then to a good degree of 
approximation, the temperature variation of optically thin dust is 
T oc 1/r*, where t — 2/(4 + n), producing the wavelength- 
dependent observed angular size 



ex 




(A6) 



Since t < \,9\ increases faster than A^, a fairly steep rise. 



© 2003 RAS, MNRAS Q0Q.ITWTT1 



Disks and Halos in PMS Stars 



11 



A3 Flux— the SED 

The flux can be obtained from eQuation lA2l bv integration over the 
observed area. At emission wavelengths, the flux at distance D is 



D2 



B,{T) m dV. 



(A7) 



Since the temperature profile of optically thin dust depends only 
on distance from the star, the dependence on the source geometry 
enters only from . 

As before, the integration is truncated by either the tempera- 
ture or the matter distribution. Whenever A > Aout at every point 
on the surface, the emission is matter bounded everywhere and 
the integration encompasses the entire source. Under this circum- 
stances Fi, oc ly^ajy, a universal SED that depends only on the 
dust properties irrespective of geometry. In particular, the spectral 
variation cr^ oc u" gives Fi, oc i/'^'^", therefore the signature of 
matter bounded emission is this SED accompanied by wavelength 
independent images; this is the expected behavior in any geom- 
etry at sufficiently long wavelengths. At A < Aout the integration 
volume is truncated by the temperature, and since the emission vol- 
ume decreases as the frequency increases, the rise of with v be- 
comes less steep than in the matter dominated regime. Therefore, 
the SED changes from F„ oc i/^"*"" at A > Aout to a flatter slope at 
A < Aout. 

The break in the slope at Aout can be used to determine the 
surface temperature Tout - Flux spectral variation shallower than i/^ 
is a clear indication of temperature-bounded emission and should 
be accompanied by an image size that increases with wavelength. 



p+t-i 



(All) 



du 



(1 +1i2){p+t)/2 

where the observed size 9\ is smaller than B when A < Aout (equa- 
tion lA6t . As long as 9 <^ dx the integration can be extended to oo 
and the brightness then decreases along any radial direction in pro- 
portion to l/6lf+*~^ 

The flux integration in eauation lA7l is similarly terminated at 
the observed boundary 9x, producing 



F.^'-^kT. 



(p+t) 



(I) 



3-(p+t) 



1 



(A 12) 



Since 9x > 9s, there are two families of SED. In the case of steep 
density distributions with p > 3 — t, the first term inside the brack- 
ets can be omitted because 3 — (p + t) < 0. Such distributions 
produce oc i/^ai, irrespective of the actual value of p. Since 
typically t ~ 0.4, this behavior applies to all cases of p > 2.6. 
On the other hand, whenever p < 3 — t the omitted term dominates 
and the SED is a broken power law. The power index switches from 
the universal 2 + n at A ^ Aout to the geometry-dependent value 
3 + n — {3 — p)/t (see also Harvey et al 1991) at A Aout- 

These results show that the SED displays a dependence on 
the density profile only when p < 2.6, and then only in the spectral 
region A Aout . All other regions of p and A produce the universal 
behavior Fi, oc v^Gu. 



A4 Spherical Geometry 

Some explicit results are easily derived in the case of spherical 
symmetry. Thanks to the scaling properties of dust radiative trans- 
fer (IE), only two properties are required to specify the geome- 
try. The first is the radial optical depth at one wavelength, say 
TV = cry J n^dr where ay is the cross-section at visual; at ev- 
ery other wavelength, Ty = q^r-v where q^ = av/aY- The second 
is the dimensionless, normalized profile of the dust density distri- 
bution 

niy) = (A8) 

n^dy 

where y = r/Rsi note that J rjdy = 1. Explicit results follow 
immediately for all power-law density profiles where 

yP i(lny) P=l 

The shell extends to the outer radius Y Rs, subtending the angu- 
lar region 9^ 9 ^ Q, where 9^ = R^/D and 6 = ¥9^. At 
scattering wavelengths 

Iv,Bca{9) = —TY L^TUvqi, \-^\ (AlO) 

"^^^^ ^ du 

(1 +m2){p+2)/2 

Whenever ^ ^ O the integration can be extended to oo, yielding 
I {9) oc l/9^^^\ the brightness decreases as a power law so long 
as the observation direction is not too close to the halo edge. At 
emission wavelengths, on the other hand. 
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